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摘要 :基于 1981 一 2013 年 华北 平原 气象 数据 , 对 华北 平原 近 30a 的 气候 时 空 变 化 趋势 与 突变 情况 进行 分 析 , 并 结合 GIMMS- 
NDVI(1981 一 2006 年 ) 和 MODIS-NDVI(2000 一 2013 年 ) 遥感 数 据 ,探讨 气 候 突 变 影 响 下 ,华北 平原 1981 一 2013 年 植被 覆盖 度 的 
空间 分 布 和 变化 特征 。 结 合生 态 学 分 区 ,从 不 同时 空 尺度 出 发 ,分 析 华 北平 原 不 同 生态 分 区 内 气候 因子 与 植被 覆盖 度 在 年 代 
际 、 年 际 和 月 变化 的 相关 关系 。 结 果 发 现 :(1)1981 一 2013 年 ,华北 平原 气温 整体 呈现 显著 上 升 趋势 (0.20"7Z10a,P<0.01) , 春 
季 气 温 的 升 高 .>15% 等 温 线 控制 范围 的 扩大 和 年 均 温 0% 等 值 线 在 华北 平原 的 消失 ,是 区 域 平均 气温 升 高 的 诱因 。 华 北平 原 
降水 整体 呈现 显著 减少 趋势 (-1.75mm/10a,P>0.05) ,其 中 秋季 降水 量 减少 过 快 ,400 一 600mm 降水 等 值 线 控制 范围 的 扩大 、 
600 一 800mm 和 800 一 1000mm 降水 等 值 线 的 范围 的 缩小 ,共同 造成 区 域 降水 量 的 减少 。 四 季 气 候 倾 向 率 的 特征 变化 敏感 区 域 
主要 位 于 北纬 35° 一 39° 之 间 。1991 一 1994 年 为 华北 平原 气候 的 突变 时 期 。(2) 华北 平原 植被 覆盖 度 总 体 呈 现 上 升 的 趋势 , 呈 
曾 加 趋势 的 面积 占 总 面积 的 55%。 人 类 活动 不 仅 加 速 了 区 域 植 被 覆盖 度 的 降低 ,也 加 剧 了 降低 速率 的 变 快 。(3) 总体 上 ,研究 
区 月 尺度 植被 覆盖 度 与 气候 因子 的 相关 性 高 于 年 尺度 的 值 。 植 被 覆盖 度 与 年 降水 量 的 偏 相 关 性 高 于 其 与 年 均 温 的 偏 相关 性 。 
年 均 温 对 农业 生态 区 和 森林 生态 区 的 植被 覆盖 度 的 影响 更 大 ,草原 生态 区 的 植被 覆盖 度 对 年 降水 量 的 依赖 性 更 强 。 在 月 尺度 
上 降水 量 对 植被 的 影响 具有 时 滞 效 应 。 气 温 对 草原 生态 区 具有 时 滞 效 应 ,降水 量 对 农业 生态 区 具有 时 滞 效 应 。(4) 华北 平原 
干 热 化 的 气候 突变 降低 了 植被 覆盖 度 的 增加 趋势 。 从 植被 覆盖 度 出 发 ,草原 生态 区 对 气候 突变 的 响应 最 明显 ;从 变化 速率 角度 
出 发 ,农田 生态 区 的 响应 最 明显 。 就 整体 而 言 ,人 类 活动 的 影响 力 还 在 持续 增强 , 且 呈 现 出 在 退化 区 的 作用 力 高 于 改善 区 的 趋 
势 。 气 候 突变 后 ,出 现 了 人 类 活动 在 植被 覆盖 度 的 改善 区 的 相对 作用 力 高 于 退化 区 的 变化 。 
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Abstract: Global climate change has resulted in significant vegetation changes in the past fifty years. The North China 
Plain, the most important grain production base of China, is experiencing prominent warming and drying. Vegetation cover, 


which is used to monitor vegetation change, is affected by climate change, predominantly temperature and precipitation. In 
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the present study, GIMMS (Glaobal Inventory Modelling and Mapping Studies) NDVI ( Normalized Difference Vegetation 
Index) ，MODIS ( Moderate-resolution Imaging Spectrometer) NDVI, and climate data, during 1981 一 2013 ，were used to 
investigate the spatial distribution and changes of vegetation. The relationship between climate and vegetation on the North 
China Plain were also analyzed on different spatial ( agricultural ，forest ，and grassland) and temporal (yearly, decadal, 
and monthly) scales. (1) It was found that the temperature increased slightly (0.20°C/10a, P < 0.01). This may be due 
to the disappearance of the 0°C isotherm, the rise of spring temperatures, and the increase of the isotherm ( >15°C ) area. 
At the same time precipitation decreased significantly (~—1.75mm/10a, P > 0.05). The expanding of the 400 一 600 mm 
precipitation contour area，the narrowing of ranges of 600 一 800 and 800 一 1000 mm ，and excessive reduction of autumn 
precipitation indicated less precipitation. The sensitive areas of characteristic climatic change were mainly located between 
the latitudes 35—39°. The climate mutation occurred during 1991 一 1994. (2) Vegetation cover increased slightly in 55% 
of the total study area, with a rate of change of 0.00039/10a. Human activities could accelerate the changes of the 
vegetation cover and affect the rate of these changes. (3) Overall, the correlation between the vegetation cover and climatic 
factors is higher for a monthly than yearly scale. The correlation analysis between vegetation cover and climate change 
showed that annual vegetation cover was positively correlated with precipitation in grassland biomes; but it showed a more 
positive correlation with temperature in the agricultural and forest biomes. In addition, vegetation cover was temporally 
sensitive to precipitation. (4) The vegetation cover increased before and after the climate variations, but the rate of increase 
slowed down. Moreover, the grassland biome had an obvious response to the climatic variations, but the agricultural biome 
showed a significant response. The effect of human activity in regions of degradation was higher than that in environmentally 
improved areas. However, after the climate abruptly changing, the effect of human activity in environmentally improvement 


areas was higher than that in the degraded regions, and the influence of human activity will continue in the future. 
Key Words: climatic change; vegetation coverage; trend analysis; residual analysis; partial correlation 


依据 IPCC 的 研究 报告 ,从 1880 到 2012 年 ,全 球 陆地 和 海洋 表面 的 平均 温度 升 高 了 0.85%C 1 。 在 过 去 的 
50a 里 ,我国 的 地 表 温 度 升 高 了 1.1%C , 升 高 的 速率 约 为 0.22/10a, 略 高 于 全 球 水 平 2 。 气 候 和 陆 气 间 相 互 作 
用 的 变化 ,通过 气温 的 升 高 ,降水 机 制 的 改变 以 及 季节 周期 性 的 波动 ,对 植物 的 生长 和 分 布 产 生 重 要 影响 5 。 
因此 ,气候 变化 引起 的 归 一 化 植被 指数 NDVI 在 一 定 范围 内 来 回 波动  。20a 以 上 的 长 期 研究 揭示 生物 和 自 
然 系统 变化 与 区 域 气温 降水 变化 相关 :5 。 因 此 区 域 气候 特征 变化 和 区 域 生 态 系统 变化 的 内 在 联系 ,以 及 两 
者 间 的 响应 机 制 成 为 了 研究 热点 。 

植被 作为 生态 系统 的 重要 组 成 部 分 ,对 诸如 大 气 系统 和 土壤 系统 等 生态 系统 之 间 的 能 量 传输 有 重要 的 影 
响 作 用 :4 。 与 植被 覆盖 度 有 正 相 关系 数 的 植被 指数 NDVI 是 指示 植被 密度 和 长 势 的 重要 指示 器 ' 。 处 在 自 
然 环境 和 人 类 活动 共同 影响 下 的 植被 覆盖 度 ,在 全 球 气 候 和 碳 循环 系统 中 扮演 重要 的 角色 。 所 以 ,植被 的 变 
化 情况 揭示 了 自然 环境 和 人 类 活动 的 相关 影响 ,以 及 两 者 之 间 的 反馈 影响 。Myneny'5 和 Tucker' ”研究 发 现 
长 时 间 序 列 的 植被 覆盖 度 呈 现 增加 的 趋势 ,特别 是 在 北半球 中 纬度 地 区 。 他 们 把 这 一 趋势 归 因 于 全 球 气候 的 
变 暖 。 一 些 研究 发 现 不 同 地 区 的 植被 覆盖 度 的 增 减 趋势 是 不 一 致 的 。Park00 和 Piao[5 研究 发 现在 欧 亚 大 陆 
增 减 趋势 是 一 致 的 。Wang'21 和 Zhang!53 研 究 发 现在 北美 洲 , 西 北部 和 东北 部 的 增 减 趋势 呈现 相反 趋势 。 已 
有 研究 都 表明 ,植被 覆盖 度 与 气候 因子 有 内 在 联系 ,尤其 是 温度 因子 和 降水 量 因子 "% 。 研 究 表明 在 中 国 的 
西北 地 区 ,植被 覆盖 度 的 变化 趋势 与 降水 量 的 变化 趋势 一 致 ,植被 指数 与 降水 量 之 间 有 显著 的 正 相 关 关系 ,与 
气温 的 相关 性 较 低 。 依 据 以 上 研究 发 现 , 植 被 覆盖 度 的 变化 呈现 明显 的 区 域 差异 性 ,与 气候 因子 关系 的 空间 
异 质 性 也 较 强 。 同 时 ,气候 因子 在 一 定 的 时 间 尺 度 内 会 有 突变 现象 发 生 号 ] 。 气 候 突 变 势 必 影 响 动 植物 生长 
和 生态 结构 变化 。 因 此 ,探讨 各 生态 系统 内 植被 覆盖 度 对 气候 变化 的 啊 应 ,对 保护 区 域 生 物 多 样 性 与 生态 安 
全 具有 重要 意义 。 
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基于 植被 上 覆盖 度 变 化 与 气候 因子 之 间 的 密切 关系 ,以 及 植被 覆盖 度 时 空 变化 的 空间 异 质 性 。 本 文选 择 华 
北平 原 为 研究 区 ,探讨 华北 平原 植被 覆盖 度 的 空间 分 布 特征 和 变化 趋势 。 并 从 气候 要 素 入 手 ,结合 生态 系统 
分 区 ,从 不 同时 空 尺 度 出 发 ,分析 了 农田 生态 区 ,草原 生态 区 和 森林 生态 区 内 植被 覆盖 度 的 年 代 际 变化 、 年 际 
变化 和 月 变化 对 降水 量 因子 和 温度 因子 的 时 空 响应 ,并 探讨 1981 一 2013 年 植被 覆盖 度 变化 对 气候 突变 的 响 
应 特征 。 


1 研究 区 地 理 概 况 及 生态 区 划分 


华北 平原 位 于 30*00' 一 40*24'N 112*48' 一 122"45 中 之 间 ,包括 河北 河南 .山东 安徽 .江苏 .北京 及 天 津 
等 5 省 2 市 的 大 部 或 部 分 地 区 ,总 面积 约 40x10'km 。 研 究 区 北 起 燕山 , 西 沿 太行 山 、 伏 牛山 , 南 抵 淮 河 干流 
及 苏 北 灌溉 总 渠 , 东 临 黄 海 和 渤海 , 除 山东 半岛 为 丘陵 外 ,大 部 分 地 区 地 势 平 坦 , 平 原 面积 约 为 33x10’km , 耕 
地 面积 约 占 总 面积 的 72%。 人 研究 区 属于 温带 大 陆 季 风气 候 区 , 光 热 资源 丰富 , 雨 热 同 季 ,降水 量 大 且 多 集中 
于 6 一 9 月。 冬小麦 . 夏 玉米 为 主要 种 植 农作物 ,种 植 制度 为 一 年 两 熟 *”。 
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1 华北 平原 地 理 位 置 和 气象 站 点 分 布 图 
Fig.l The distribution of the meteorological station 
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参考 Wu 等 中 的 研究 成 果 , 将 研究 区 分 为 8 个 地 
理 分 区 (图 2) 。 从 南 到 北 按照 编码 顺序 依次 为 :淮南 与 
长 江 中 下 游 自 然 区 (1) 汉中 盆地 自然 区 (2) .华北 平原 
自然 区 (3) 、 鲁 中 山地 丘陵 自然 区 (4) .辽东 胶东 山地 丘 
陵 区 (5) ,华北 山地 丘陵 自然 区 (6) .大 兴安 岭南 部 自然 
区 (7) 和 内 蒙古 高 原平 原 东部 自然 区 (8)。 为 探讨 不 同 
生态 系统 分 区 植被 覆盖 度 的 变化 特征 ,将 已 有 自然 分 区 
按照 生态 系统 类 型 合并 为 3 个 大 类 分 区 ;森林 生态 系统 
区 、 草 原生 态 系 统 区 和 农业 生态 系统 区 。 结 果 如 图 2 所 
示 。 草 原生 态 区 主要 分 布 于 河北 省 的 西北 部 ,包括 张 北 
地 区 和 木兰 围场 。 该 生态 区 地 处 我 国 温带 半 湿润 地 区 
与 半 干 旱地 区 的 交织 带 , 主要 以 温 性 草 多 草原 .山地 草 
甸 和 温 性 草原 为 主 2] 。 森 林 生 态 区 包括 燕山 山脉 . 太 
行 山 - 伏 牛山 山脉 和 山东 半岛 的 丘陵 山地 ,植被 类 型 从 图 2 研究 区 生态 区 分 区 图 
山 枇 至 山顶 ,大 体 分 布 落叶 林 、 冰 叶 针 叶 混交 林 、 针 叶 林 Fig.2 The distribution of the ecological area 
和 高 山 灌木 草 从 ,主要 以 松林 、 侧 柏林 、 栎 林 为 主 [2。 
农田 生态 区 分 布 区 域 为 除了 草原 生态 区 和 森林 生态 区 以 外 的 区 域 。 华 北平 原 的 农业 以 一 年 两 熟 为 主 ,农作物 
中 冬小麦 .玉米 占 绝 对 优势 。 


2 数据 准备 与 研究 方法 


2.1 数据 准备 与 一 致 性 检验 

美国 马里 兰 大 学 提供 的 CGIMMS-NDVI 数据 集 包 含 1981 一 2006 年 半月 最 大 合成 数据 ,空间 分 辨 率 为 8kmx 
8km, 数 据 经 过 辐射 校正 和 几何 校正 .大 气 校正 、 除 云 等 处 理 ,GIMMS 数据 集 的 质量 优 于 其 他 NDVI 数据 
集 ' 引 。GIMMS 数据 集 是 目前 最 长 时 间 序 列 的 NDVI 数据 ,与 其 他 高 分 辩 率 数据 集 之 间 具有 较 好 的 相关 性 ,已 
广泛 应 用 在 全 球 和 区 域 尺度 植被 动力 学 研究 中 ' 关 5 。NASA 提供 的 MODIS-NDVI 数据 集 包含 2000 一 2013 年 
华北 平原 MOD13Q1 级 产品 ,时 间 分 辩 率 为 16d, 空 间 分 辩 率 为 250mx250m。 应 用 MODIS Reprojection Tools 
(MRT) 软 件 对 数据 进行 投影 转换 .数据 拼接 等 处 理 。 

在 获取 植被 覆盖 度 的 方法 中 ,最 大 值 合 成 法 MVC ( Maximum Value Composites ) 使 用 较为 广泛 , 它 是 基于 
NDVI 名 数 据 合成 获得 月 度 和 年 度 NDVI 数据 的 方法 。 具 体 公式 如 下 051 : 

NDVL = Max( NDVI,) 

式 中 ,NDVI, 指 第 i 个 月 或 者 第 i 年 的 NDVI，NDVI; 指 第 i 月 的 第 j 个 15 天 的 NDVI 数据 或 第 i 年 的 第 j 月 的 
NDVI 数据 。 

由 于 GIMMS 和 MODIS 数据 采用 不 同 的 传感器 ,因此 需要 对 两 者 进行 一 致 性 检验 ”3。 本 文采 用 的 
CIMMS-NDVI 数据 的 时 间 序 列 为 1981 一 2006, MODIS-NDVI 数据 的 时 间 序 列 为 2000 一 2013 年 ,依据 2000 一 
2006 年 的 重 车 数 据 ,对 年 最 大 NDVI 数据 进行 相关 性 分 析 , 得 到 两 者 相关 系数 为 0.801 ,在 0.05 置信 水 平 上 显 
著 相 关 ,说 明 两 种 数据 的 年 度数 据 在 全 区 域 尺度 上 具有 显著 一 致 性 ,二 者 的 线性 回归 方程 为 :NDVIcwws = 
0.5886xNDVIvwoms+0.0965(*=0.6417,n=7,P<0.05)。 采 用 同样 方法 对 两 者 的 月 尺度 数据 进行 一 致 性 检验 ， 
两 者 月 尺度 数据 的 相关 系数 为 0.957, 在 0.01 置信 水 平 上 通过 检验 ,线性 回归 方程 为 :NDVIcwws = 0.781 
NDVIwowms 一 0.0469(r*=0.9161,n=7,P <0.01) 。 通 过 以 上 回归 方程 ,结合 2007 一 2013 年 MODIS-NDVI 数据 插 
补 2007 一 2013 年 度 和 月 度 GIMMS-NDVI 数据 ,将 实验 数据 的 时 间 跨 度 延 伸 为 1981 一 2013 年 。 
2.2 ”突变 检测 

气候 突变 是 普遍 存在 于 气候 系统 中 的 重要 现象 。 气 候 突变 的 普 适 性 定义 :气候 从 一 种 稳定 态 ( 或 稳定 持 
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续 的 变化 趋势 ) 跳跃 式 地 转变 到 另 一 种 稳定 态 (或 稳定 持续 的 变化 趋势 ) 的 现象 ,表现 为 气候 在 时 空 上 从 一 个 
统计 特性 到 另 一 个 统计 特性 的 急剧 变化 。 在 气候 突变 的 检测 方法 中 ,Mann-Kendall 非 参数 检测 法 的 检测 范围 
宽 、 人 为 性 少 ,定量 化 程度 高 ,不 需要 样本 遵从 一 定 的 分 布 ,也 不 受 少数 异常 值 的 干扰 |。 

Mann-Kendall 检测 的 基本 原理 如 下 :气候 序列 设 为 平稳 序列 (i=1,2,3,…,n) ,表示 第 i 个 样本 (1<j<i) 
的 累计 数 ,定义 统计 量 : 


di = py Mm; 
在 原 序列 的 随机 独立 等 的 假设 下 ,均值 方差 分 别 为 : 
var(di;)=k(k— 1)(2k + 5)/72 
di 标准 化 后 ， 
UP = (di —- E(di))/ Vvar(di) 
满足 标准 正 态 分 布 ,其 概率 a=prob(1U1>1UF, ) 可 以 通过 计算 或 者 查 表 获 得 。 给 出 显著 性 水 平 a( 取 a= 
0.05,=1.96) , 若 1VFI>U, 则 表明 该 气候 序列 呈现 显著 的 变化 趋势 。 
将 气候 序列 (i =1,2,3,…,n) 的 逆序 列 ,按照 以 上 过 程 计 算 一 次 。 同 时 使 1UB,1=-UF, ,k=nn-1,...,1， 
VB,=0。 若 1UFI>0, 此 时 序列 旦 上升 趋 势 ; 若 1UFI<0 序列 呈 下 降 趋 势 。 通 过 信 度 检验 可 知 其 是 否 有 变化 趋 
势 。 当 曲线 超过 置信 水 平 线 时 ,表明 序列 的 上 升 (下 降 ) 趋 势 显著 . 若 UF 和 UB 两 条 曲线 出 现 交 点 , 且 交 点 位 
于 置信 区 间 内 ,该 交点 对 应 的 时 刻 即 序列 突变 开始 的 时 刻 。 
2.3 植被 覆盖 度 变化 趋势 
基于 像 元 的 植被 覆盖 度 的 变化 趋势 ,可 以 通过 一 元 线性 回归 分 析 得 出 ,计算 公式 为 5 : 


nx Si- 
式 中 ,0,,. 为 变化 趋势 的 斜率 ,n 为 监测 的 年 数 , C, 为 第 ; 年 的 年 最 大 植被 覆盖 度 。 植 被 覆盖 年 际 间 变化 的 显 
著 性 可 以 通过 年 时 间 序 列 和 植被 覆盖 度 的 相关 关系 获得 。 正 值 代表 植被 覆盖 度 上 升 , 负 值 代表 植被 覆盖 度 下 
降 。 变 化 趋势 的 显著 性 检验 采用 下 检验 。 统 计量 计算 公式 为 : 
N-2 
0 

其 中 , U = 之 ($, - 7) ”为 误差 平方 和 , 0 = 之 (y; -分 ) 为 回归 平方 和 ,y, 为 第 i 年 覆盖 度 真实 观测 值 ,为 
覆盖 度 的 10 年 平均 值 。 根 据 检验 结果 将 趋势 变化 分 为 5 个 等 级 . 极 显著 减少 9,,,.<0,P<0.01); 显著 减少 
bs<0,0.01<P<0.05) ;变化 不 显著 (P>0.05); 显著 增加 (0,,,.>0,0.01<P<0.05); 极 显 著 增 加 9,,.>0,P< 
0.01) 。 
2.4 气候 因子 与 植被 覆盖 度 的 相关 关系 

地 理 系统 是 一 种 多 要 素 的 复杂 巨 系 统 ,其 中 一 个 要 素 的 变化 必然 影响 到 其 他 要 素 的 变化 。 在 多 要 素 系统 
中 ,研究 一 个 要 素 对 另外 一 个 要 素 的 相关 性 时 , 暂 不 考虑 其 他 要 素 的 影响 , 则 称 为 偏 相关 '”i。 偏 相关 系数 的 
计算 公式 为 : 


F=Ux 


(1 -7.)(1 -7.) 
式 中 , 7, 为 变量 z 固定 后 变量 x 与 变量 y 的 偏 相关 系数 ; r, 为 变量 * 与 变量 y 的 相关 系数 ; r, 为 变量 y 与 变 
量 z 的 相关 系数 ; +, 为 变量 x 与 变量 z 的 相关 系数 。 
其 中 ,zy 和 z 之 间 的 相关 系数 计算 公式 为 : 


rx = 
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式 中 ,7 为 x 与 ;之 间 的 相关 系数 ,x wy 分 别 为 第 i 年 第 j 月 的 研究 变量 值 ;x, wy 分 别 为 变量 在 各 时 间 尺 度 的 
均值 。 
偏 相关 系数 的 显著 性 检验 ,一 般 采 用 ,检验 法 。 其 统计 量 计算 公式 为 : 


wV7 一 1 一] 


M1- Ta 
式 中 , ra 为 偏 相 关系 数 ;n 为 样本 数 ;m 为 自 变 量 个 数 。 显 著 水 平 的 临界 值 1, 可 以 通过 查询 ;分布 表 获 
得 , 若 i>i,, 则 表明 偏 相关 显著 ;车 i<i, , 则 偏 相关 不 显著 。 
2.5 ”和 残 差分 析 

气候 变化 对 植被 覆盖 度 变化 产生 重要 影响 ,降水 因子 和 和 气温 因子 是 气候 变化 影响 中 的 重要 因子 。 通 过 降 
水 因子 和 气温 因子 拟 合 出 NDVI 的 预测 值 , 即 为 气候 变化 的 影响 量 。 计 算 其 与 NDVI 真实 值 之 间 的 差 ,可 作为 
人 为 因素 的 影响 ,该 方法 即 为 残 差 分 析 , 已 经 得 到 广泛 应 用 [3 。 

2=NDVI.-NDVL,。 

式 中 ,NDVIL,, 为 NDVI 真实 值 , NDVI 为 NDVI 预测 值 ,s>0 ,说 明 人 类 活动 产生 正面 影响 ;se<0, 说 明 人 类 活动 
产生 负面 影响 ;se=0, 说 明 人 类 活动 影响 比较 微弱 。 


一 


3 ”结果 与 讨论 


3.1 华北 平原 气候 变化 特征 与 突变 分 析 

1981 一 2013 年 华北 平原 气温 和 降水 量 的 年 际 和 四 季 变 化 趋势 如 网 3 所 示 , 表 1 为 统计 的 气温 和 降水 量 
的 年 际 和 四 季 的 气候 倾向 率 。 从 图 3 中 可 见 ,华北 平原 年 均 温 呈现 整体 上 升 的 趋势 ,其 气候 倾向 率 约 为 
0.26%7Z10a(P<0.01) 。 和 气温 的 季节 性 变化 整体 呈现 升 高 趋势 ,只 是 幅度 略 有 不 同 ,春季 (0.34%7Z10a ) 增 温 最 
快 ,冬季 (0.09%7Z10a) 最 慢 。 春 季 均 温 的 气候 倾向 率 大 于 年 均 温 ,说 明 春 季 增 温 对 年 均 温 的 增长 贡献 大 于 其 
他 季节 。 该 研究 结果 与 赵 少 华 的 研究 结果 一 致 :1 。 降 水 量 呈 现 减少 的 趋势 ,气候 倾向 率 约 为 -18.13mm/10a 
(P>0.05)。 与 气温 的 季节 性 整体 升 高 趋势 不 同 ,降水 量 季 节 性 差异 更 加 明显 。 降 水 量 在 秋季 减少 最 快 
(-10.70 mm/10a) ,夏季 减少 最 慢 (-0.35 mm/10a) ,冬季 呈现 略微 增长 的 趋势 ,气候 倾向 率 为 1.4mm/10a。 虽 
然 冬季 的 降水 呈 增 加 趋势 ,但 增加 幅度 小 , 远 不 及 其 他 季节 降水 量 的 减少 幅度 。 从 气温 与 降水 量 的 变化 趋势 
可 见 ,华北 平原 的 气候 变化 有 向 干 热 化 发 展 的 趋势 。 


表 1 1981 一 2013 年 气候 因子 气候 倾向 率 统计 结果 


Table 1 The climate tendency rate statistical results of climate factor from 1981 to 2013 


时 间 尺 度 气温 因子 Temperature 降水 量 因 子 Precipitation 

Time scale 一 元 线性 方程 Linear equation R 一 元 线性 方程 Linear equation R 
春 Spring y=0.03x+12.89 0.39 y=-0.84x+167.96 0.18 
夏 Summer y=0.02x+24.87 0.18 y=-0.03x+433.95 0.05 
秋 Autumn y=0.018x+13.73 0.22 y=—1.07x+164.32 0.23 
冬 Winter y=0.009x+0.22 0.08 y=0.14x+54.986 0.18 
年 Year y=0.026x+12.677 0.46 y=—1.80x+815.55 0.18 


图 4 为 华北 平原 年 降水 量 和 年 均 温 的 气候 倾向 率 的 空间 分 布 图 。1981 一 2013 年 间 ,华北 平原 82 个 气象 
站 点 的 年 均 温 气 候 倾 向 率 大 部 分 呈现 增 温 变化 趋势 , 增 温 趋势 为 0.32%/10a。 增 温 趋势 较 高 区 域 位 于 河北 省 
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图 3 华北 平原 气温 与 降水 量 数据 的 变化 趋势 


Fig.3 The change trend of temperature and precipitation 


西部 太行 山地 区 以 及 山东 省 的 威海 (0.64%C/10a) 和 青岛 (0.67%CA10a)。 河 北 承 德 山 东 东 营 站 点 年 均 温 的 气 
候 倾向 率 呈 现 降低 趋势 ,气候 倾向 率 分 别 为 -0.06%Z10a 和 -0.28%7Z10a。 河 北 省 的 气候 倾向 率 均 值 高 于 山东 
省 ,山东 省 的 值 高 于 河南 省 。 说 明 在 全 球 增 暧 背景 下 ,华北 平原 气温 的 空间 变化 具有 一 定 的 纬度 地 带 性 , 即 : 
气温 的 气候 倾向 率 随 纬度 升 高 而 升 高 。 春 季 夏季 和 秋季 的 气候 倾向 率 分 布 特征 从 南 到 北 呈 现 高 一 低 一 高 的 
趋势 , 低 值 区 域 分 布 也 相似 ,主要 呈 带 状 沿 河南 .山东 和 河北 三 省 交界 地 区 向 东 直 到 渤海 一 带 分 布 。 冬 季 从 南 
到 北 呈 现 逐 渐 增 大 的 趋势 。 华 北平 原 年 降水 量 的 气候 倾向 率 为 -1.75%CZ10a ,整体 呈 减 少 趋势 ,减少 幅度 最 快 
的 站 点 为 东营 。 年 降水 量 增 加 的 区 域 主 要 位 于 南阳 一 郑州 一 菏泽 一 线 的 南 侧 。 春 、 秋 两 季 呈 现 南部 减少 , 北 
部 增加 的 趋势 , 夏 、 冬 两 季 呈 现 南 部 增加 北部 减少 的 趋势 。 季 节气 候 倾 向 率 增 减 特征 变化 敏感 区 域 主要 位 于 
35° 一 39°N 之 间 。 

分 别提 取 各 年 代 际 的 气候 因子 的 面积 要 素 , 探寻 气温 和 降水 量 因子 各 量 级 在 不 同年 代 际 的 空间 分 布 变 
化 。 气 温 因子 按照 <0 .0 一 5 .5 一 10 10 一 15% 、>15%C 五 类 统计 ,降水 量 因子 按照 <400 .400 一 600 .600 一 800、 
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4 华北 平原 年 降水 量 气候 倾向 率 (a) 与 年 均 温 气候 倾向 率 (b) 空间 分 布 


Fig.4 The spatial distribution of climate tendency rate of precipitation (a) and temperature (b) 


800 一 1000 .1000 一 1500 、>1500mm 六 类 统计 。 

从 图 5 中 可 以 发 现 ,1960 一 2010 年 ,华北 平原 大 部 分 地 区 都 被 10 一 15% 范围 的 温度 控制 (多 年 均值 面积 
占 比 78.52% ) ,其 次 是 5 一 10%C 的 控制 范围 (15.49%) ,>15%C 的 控制 范围 排 在 第 三 位 (5.08% ) ,<0% 和 0 一 5% 
的 平均 面积 占 比 均 未 超过 1%。 各 量 级 区 间 面 积 占 比 在 不 同年 代 际 的 变化 各 不 相同 。1980s 是 各 温度 范围 控 
制 面积 变化 的 特征 年 代 ,>15% 的 面积 范围 在 1990s 超过 了 5 一 10%C 的 面积 范围 。 总 体 上 ,>15% 温度 的 控制 
范围 在 变 大 ,其它 等 温 线 控制 的 范围 ,从 1960s 到 2010s 面积 都 在 减 小 。 结 合 研究 区 年 均 温 的 整体 升 高 趋势 ， 
可 以 得 出 结论 :一 方面 ,全 球 变 暖 背景 下 的 华北 平原 各 站 点 的 年 均 温 大 部 分 呈 升 高 趋势 ; 另 一 方面 ,>15%C 等 
温 线 控制 范围 的 扩张 ,也 是 年 均 温 升 高 的 因素 之 一 。 从 图 5 可 见 ,华北 平原 降水 量 400 一 600mm 等 值 线 ( 面积 
占 比 33.97% ) 和 600 一 800mm (面积 占 比 35.06%) 等 值 线 控制 范围 的 平均 面积 相当 ,其 次 是 降水 量 1000 一 
1500mm 等 值 线 控制 范围 (面积 占 比 20.15%)。<400mm 等 值 线 控制 面积 和 1000 一 1500mm 的 控制 面积 的 平 
均值 相当 ,分 别 为 4.51% 和 6.87%。 研 究 区 降水 量 >1500mm 的 面积 占 比 多 年 来 一 直 低 于 1% ,但 也 呈现 略微 
下 降 趋 势 。 总 体 上 ,全 区 大 部 分 地 区 的 降水 量 范 围 为 400 一 600mm 和 600 一 800mm , 两 种 类 型 斑 块 面积 平均 占 
比 之 和 达到 79%。400 一 600mm 等 值 线 控制 的 面积 整体 呈现 增加 趋势 ,其 余 等 值 线 的 控制 范围 都 呈现 下 降 或 
者 维持 不 变 的 态势 。1970s 一 1980s 之 间 ,400 一 600mm 降水 等 值 线 控制 范围 超过 了 600 一 800mm 控制 范围 ,成 
为 最 大 面积 斑 块 。1980s 是 各 降水 量 等 值 线 控制 面积 变化 的 特征 年 代 。 

华北 平原 的 气候 环境 在 过 去 的 30 多 年 内 经 历 了 突变 ,尤其 是 气温 的 突变 比较 明显 (图 6) 。M-K 检测 结 
果 显 示 :1981 一 2013 年 间 ,UF 和 UB 曲线 在 1994 年 相交 , 且 交 点 在 置信 区 间 内 。 累 积 距 平 统计 结果 呈 “V? 字 
形 ,转折 点 出 现在 1993 年 ,表明 华北 平原 气温 经 历 了 从 1981 一 1993 年 之 间 的 降温 过 程 和 1994 一 2009 年 之 间 
的 增 温 过 程 ,2009 一 2013 年 又 有 小 幅 回 落 。 结 合 两 种 检测 方法 结果 分 析 得 出 :1981 一 2013 年 ,华北 平原 气温 
由 冷 到 暧 的 转变 时 间 出 现在 1993 一 1994 年 间 。 华 北平 原 年 降水 量 的 M-K 检测 结果 显示 (图 7) ,UF 和 UB 曲 
线 在 置信 水 平 内 有 5 个 交点 。 依 据 累 积 距 平 的 统计 结果 分 析 ,1982 一 1991 年 间 ,降水 量 累积 距 平 呈 现 增多 的 
变化 过 程 ;1991 一 2013 年 间 ,累积 距 平 变 为 减少 趋势 ,降水 量 累 积 距 平 的 最 大 值 出 现在 1991 年 。 结 合 两 种 检 
测 结果 分 析 得 出 :在 1981 一 2013 年 间 ,华北 平原 降水 由 多 到 少 的 转变 时 间 发 生 在 1991 年 前 后 。 结 合 年 均 温 
和 降水 量 的 突变 研究 结果 得 出 :1991 一 1994 年 间 是 气候 因子 活跃 变化 时 期 ,该 时 期 为 突变 特征 时 期 。 
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图 5 华北 平原 气温 与 降水 量 各 统计 量变 化 趋势 图 


Fig.S The change trend chart of temperature and the precipitation in Huabei Plain 
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图 6 年 均 温 M-K 检测 与 累积 距 平 图 


Fig.6 M-K test and Cumulative anomaly of temperature 


3.2 ”植被 覆盖 度 空间 分 布 与 变化 趋势 特征 
图 8 是 华北 平原 1981 一 2013 年 植被 覆盖 度 空 间 分 布 格局 。 人 研究 区 年 平均 植被 覆盖 度 为 0.6471。 城 市 
群 、 较 大 的 水 面 和 海岸 线 地 区 的 植被 覆盖 度 相 对 稀 玻 ,甚至 没有 植被 覆盖 。 植 被 覆盖 度 较 低 和 中 等 水 平 植被 
善 度 区 域 , 主 要 分 布 在 河北 省 西北 部 半 干 旱 的 农 牧 交 织带 以 及 河南 省 和 山东 省 的 低 山 丘陵 地 区 。 植 被 覆盖 
度 较 高 的 区 域 主要 集中 在 华北 平原 中 部 ,包括 河北 平原 .河南 省 的 中 东部 .山东 省 西部 以 及 安徽 北部 地 区 。 植 
被 覆盖 度 最 高 区 域 集中 分 布 在 燕山 和 秦岭 东 延 的 伏 牛 山 山 脉 中 的 几乎 不 受 人 类 影响 的 区 域 。 草 原生 态 区 多 
年 平均 植被 覆盖 度 约 为 0.5504。 和 森林 生态 区 平均 植被 覆盖 度 约 为 0.6526。 农 田 生 态 区 平均 植被 覆盖 度 约 为 
0.6520 ,水 平 与 和 森林 生态 区 差别 不 大 。 
从 月 尺度 植被 覆盖 度 变 化 趋势 图 (图 9) 中 可 以 看 出 ,在 植被 生长 季 , 草 原生 态 区 和 森林 生态 区 的 植被 覆 
盖 度 在 植物 生长 期 内 表现 出 先 增 加 后 下 降 的 趋势 。 两 者 出 现 峰值 均 在 8 月 份 ,分 别 为 0.541 和 0.549。 草 原 
生态 区 的 植被 覆盖 度 整体 低 于 森林 生态 区 。 农 业 生 态 区 的 植被 覆盖 度 在 植物 生长 季 的 变化 相对 复杂 ,6 月 份 
植被 覆盖 度 出 现 明 显 的 降低 ,随后 持续 增高 ,在 8 月 出 现 峰值 0.59。 华 北平 原 农作物 以 冬小麦 为 主 , 6 月 是 小 
麦收 获 的 季 , 大 面积 裸露 的 农田 造成 了 该 生态 区 植被 覆盖 度 的 突然 降低 。4 一 5 月 的 峰值 期 ,农业 生态 区 的 植 
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图 7 年 均 降水 量 M-K 检测 与 累积 距 平 图 


Fig.7 M-K test and Cumulative anomaly of precipitation 


被 类 型 主要 是 以 小 麦 为 主 。8 月 的 峰值 期 主要 以 玉米 高梁 以 及 经 济 作物 为 主 , 且 生态 区 内 栽培 果树 等 都 达 
到 生长 季 峰 值 。 植 被 类 型 的 不 同 ,是 造成 8 月 份 峰值 高 于 4 一 5 月 峰值 的 主要 原因 。 
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图 8 植被 覆盖 度 空间 分 布 图 9 植被 生长 期 内 各 生态 区 覆盖 度 变 化 趋势 
Fig.8 Spatial patterns of vegetation coverage Fig.9 Variation trend of coverage of ecotope 


图 10 为 基于 像 元 的 植被 覆盖 度 变化 趋势 空间 分 布 图 。 华 北平 原 植 被 覆 善 度 呈 现 增 加 趋势 ,增加 速率 为 
0.00039/10a ,呈现 增加 趋势 的 区 域 占 总 面积 的 55% , 主要 位 于 河北 省 西北 部 山区 .山东 省 西南 部 河南 省 东部 
以 及 安徽 省 北部 地 区 。 呈 现 减少 趋势 的 区 域 主要 集中 于 人 类 活动 相对 频繁 的 地 区 ,包括 河南 省 中 部 和 西部 ， 
山东 省 北部 ,以 及 人 口 密集 京 津 唐 地 区 。 在 各 生态 区 的 趋势 变化 中 ,森林 生态 区 的 增加 趋势 最 快 , 约 为 
0.000034/10a; 农 业 生 态 区 的 增加 速率 最 慢 , 约 为 0.00024/10a, 草 原生 态 区 增加 的 速率 约 为 0.0029/10a。 

3.3 气候 因子 与 植被 覆盖 度 的 关系 
3.3.1 不 同年 代 际 的 水 热 条 件 对 植被 覆盖 度 的 影响 

植被 生长 受到 气象 要 素 .土地 利用 和 C0, 的 等 因素 的 多 重 影响 ,其 中 气温 和 降水 是 植被 生长 的 关键 性 因 
素 !? ,水 热 条 件 是 决定 植被 空间 分 布 及 其 变化 的 主要 非 生物 因素 福 ] 。 同 时 , 地表 植被 覆盖 度 的 变化 也 影响 
局 部 气候 。 为 探讨 不 同 水 热 组 合 与 植被 覆盖 度 之 间 的 相互 作用 关系 ,选择 1980s .1990s 和 2010s 3 个 时 期 
的 相关 数据 ,分别 计算 温度 累积 距 平 .降水 量 距 平 百分率 和 植被 覆盖 度 的 变化 趋势 ,探寻 不 同 水 热 组 合影 响 下 
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的 植被 覆盖 度 变 化 情况 。1980s .1990s 和 2010s 3 个 时 期 的 气温 和 降水 相 较 于 多 年 均值 的 锰 亏 情况 如 表 2 所 
示 :1980s 年 处 于 气温 偏 低 、 降 雨 偏 多 的 水 热 状况 ,1990s 处 于 气温 偏 高 ,降雨 偏 多 的 水 热 状 况 ,2010s 处 于 气温 
偏 高 降雨 偏 少 的 水 热 状况 。 


0.25 『 
气温 
3 
己 0.15 下 
b= 
加 本 
捍 
六 
SS 
名 写 0.05 上 
和 抠 w 
: 入 多 ky 
司 型 8 
三 中 -0.05 上 
SS 5 9 
nh 所 
[aad 
-0.15 上 
0 100 200km 
| | 
115°00'00” 120°00'00"E 1 | 1 | | 1 
变化 趋势 0 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 
1 
国 -0.011 一 -0.004 Ey] -0.001 一 0.001 Ml 0.004 一 0.018 年 份 Y 
国 -0.004 一 -0.001 Bl 0.001 一 0.004 
10 植被 覆盖 度 变 化 趋势 空间 分 布 图 11 1990s 研究 区 气候 因子 偏差 分 析 
Fig.10 Dynamic of vegetation coverage Fig.11 Deviation analysis of climatic factors in 1990s 
表 2 不 同年 代 际 气温 和 降水 的 增 减 趋势 
Table 2 Trends of temperature and precipitation over decadal 
时 间 尺 度 温度 距 平 百分比 Temperatur departure percentage 降水 量 距 平 百分比 Precipitation departure percentage 
Time scale 1980s 1990s 2010s 1980s 1990s 2010s 
春季 Spring —0.5034 —0.151 0.791 —0.02146 0.101 —0.119 
夏季 Summer -0.415 0.072 0.371 -0.015 0.014 -0.003 
秋季 Autumn —0.348 0.083 0.279 0.140 -0.105 -0.010 
冬季 Winter —0.676 0.407 0.191 -0.179 -0.010 0.219 
年 Year -0.512 0.096 0.449 0.0002 0.008 -0.011 


在 1980s, 华 北平 原 植被 覆盖 度 呈 现 整 体 增加 的 趋势 ,速率 约 为 0.03/10a。 农 田 生 态 区 (0.08/10a) 和 森林 
生态 区 (0.04/10a) 呈 增加 趋势 ,草原 生态 区 呈现 减少 趋势 ( -0.03/10a) 。 该 年 代 际 的 温 低 水 多 的 水 热 条 件 对 
整体 以 及 农田 和 森林 生态 区 的 植被 覆盖 度 呈 现 增加 的 影响 ,对 草原 生态 区 则 呈现 减少 的 影响 。 

在 1990s ,华北 平原 植被 覆盖 度 的 呈现 整体 减少 趋势 ,速率 约 为 -0.02/10a。 农 业 生态 区 ( -0.05/10a) 和 和 森 
林 生 态 区 ( -0.03/10a ) 呈 减 少 趋势 ,农业 生态 区 的 减少 速率 最 快 ,草原 生态 区 呈现 增加 的 趋势 (0.02/10a) 。 结 
合 该 时 期 的 年 均 温和 年 降水 量 的 距 平 百分比 可 以 发 现 , 雨 热 不 同步 抑制 了 植被 的 生长 (图 11) 。 草 原生 态 区 
的 植被 覆盖 度 呈 现 增 加 趋势 , 主要 得 益 于 5% 等温 线 和 400mm 降水 等 值 线 的 北 移 。 以 400mm 降水 等 值 线 为 
例 ,多 年 均 位 于 草原 生态 区 的 边缘 ,而 该 时 段 却 几乎 覆盖 了 整个 草原 生态 区 。 相 较 于 整个 研究 区 的 雨 热 不 均 ， 
平原 草原 区 降水 的 增多 结合 适宜 的 温度 ,导致 植被 覆盖 度 呈 现 增加 趋势 。 

在 2010s ,华北 平原 植被 覆盖 度 呈 现 整体 增加 趋势 ,速率 约 为 0.02/10a。 农 田 生 态 区 (0.02/10a) 、 和 森林 生 
态 区 (0.07/10a) 和 草原 生态 区 (0.13/10a) 均 呈现 增加 趋势 。 研 究 区 的 干 热 化 趋势 已 经 呈现 ,区 域内 植被 覆盖 
度 呈 现 整体 增加 的 趋势 。 虽 然 降 水 量 整体 有 所 降低 ,但 是 与 多 年 均值 的 空间 分 布 相 比 ,没有 呈现 明显 差异 。 
说 明 区 域内 各 生态 区 的 气温 和 降水 量 达 到 了 较 好 的 水 热 组 合 条 件 ， 
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综合 1981 一 2010 的 华北 平原 不 同年 代 际 的 植被 覆盖 度 对 不 同 水 热 组 合 的 分 析 结 果 得 出 ;在 基于 年 代 际 
的 时 间 尺 度 研 究 中 ,气温 的 增加 ,对 华北 平原 农业 生态 区 的 影响 体现 为 植被 覆盖 度 的 稳定 增加 ,但 是 速率 在 不 
断 降 低 。 降 水 量 的 减少 对 农业 生态 区 植被 覆盖 度 的 影响 相对 较 小 。 草 原生 态 区 植被 覆盖 度 的 变化 趋势 在 一 
定 程 度 上 与 气温 的 变化 趋势 一 致 。 当 气温 偏 高 降雨 偏 多 时 ,森林 生态 区 的 植被 覆盖 度 呈 减少 趋势 ,其 余 的 水 
热 组 合 方式 植被 覆盖 度 均 呈现 增加 趋势 。 
3.3.2 ”气候 突变 前 后 植被 覆盖 度 变 化 特征 

在 3.1 节 得 出 结论 1991 一 1994 年 为 研究 时 段 华北 平原 气候 因子 突变 的 特征 时 段 。 因 要 对 气候 突变 前 后 
植被 覆盖 度 对 气候 因子 响应 情况 进行 比 对 ,因此 暂 将 华北 平原 的 气候 突变 年 选 在 1994 年 ,此 时 降水 量 因子 和 
气温 因子 均 已 经 历 了 突变 变化 。 

将 1981 一 1994 年 和 1995 一 2013 年 两 个 时 间 段 的 遥感 影像 分 别 计算 ,得 出 基于 像 元 的 植被 覆盖 度 分 布 和 
变化 速率 统计 结果 图 (图 12) 。 华 北平 原 年 均 植 被 覆盖 度 在 气候 突变 前 后 分 别 为 0.662 和 0.663。 结 合 各 生 
态 区 年 均 植 被 覆盖 度 的 统计 结果 图 发 现 ,草原 生态 系统 区 在 气候 突变 后 由 0.544 增加 到 0.590 , 显著 增幅 约 为 
8%。 农 田 生 态 系 统 区 由 0.668 增加 到 0.669。 和 森林 生态 系统 区 由 0.668 减少 到 0.663。 说 明 植 被 覆盖 度 在 研 
究 区 内 部 ,空间 异 质 性 较 强 ,气候 突变 对 草原 生态 区 的 植被 覆盖 度 的 影响 比 其 他 生态 区 显著 。 气 候 突 变 前 后 
植被 覆盖 度 的 平均 增加 趋势 为 0.018/10a 和 0.006/10a, 增 加 幅度 降低 了 0.012/10a。 在 各 生态 区 中 ,降低 幅度 
一 ”最 快 的 是 农田 生态 区 (降低 了 0.016/10a) ,其 次 是 森林 生态 区 (降低 了 0.005/10a) ,草原 生态 区 变化 最 小 ( 降 
号 低 了 0.001/10a) 。 在 气候 突变 前 ,增加 速率 最 快 的 是 农田 生态 区 ,气候 突变 后 草原 生态 区 植被 覆盖 度 增加 速 
入 ” 率 最 快 。 森 林 生 态 区 的 速率 相对 较 慢 。 气 候 突变 对 农田 生态 区 的 植被 覆盖 度 变 化 速率 影响 最 大 ,其 次 是 森林 
OO 生态 区 ,草原 生态 区 几乎 不 受气 候 突变 的 影响 。 可 见 , 从 植被 覆盖 度 高 低 的 角度 出 发 ,气候 突变 对 草原 生态 区 
一 影响 最 大 ;从 变化 速率 角度 出 发 ,对 农田 生态 区 的 影响 最 大 。 
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12 气候 突变 前 后 各 生态 区 植被 覆盖 度 变 化 情况 统计 结果 


Fig.12 The results of climate change effects on Vegetation coverage of ecological zones 


3.3.3 植被 覆盖 度 与 气候 因子 的 相关 性 分 析 

图 13 为 年 际 植被 覆盖 度 与 年 均 温 和 年 降水 量 的 偏 相关 性 空间 分 布 结果 。 如 图 显示 :在 1981 一 2013 年 之 
间 ,华北 平原 年 际 植被 覆盖 度 与 年 降水 量 的 整体 平均 偏 相 关系 数 为 0.018 ,年 植被 覆盖 度 与 年 温度 的 整体 平均 
偏 相 关系 数 为 0.006 ,其 中 呈正 相关 的 像 元 分 别 占 总 像 元 面积 的 59% 和 51% 。 总 体 来 说 ,植被 覆盖 度 与 降水 
量 的 偏 相关 性 ,高 于 其 与 年 均 温 的 偏 相关 性 。 两 者 在 空间 分 布 上 都 表现 出 明显 的 空间 异 质 性 。 

农业 生态 区 的 植被 覆盖 度 与 年 均 温 和 年 降水 量 的 平均 相关 系数 为 0.01 和 -0.001 ,与 之 对 应 的 呈正 相关 
的 区 域 面积 占 分 区 面积 的 52% 和 48% ,表明 在 农业 生态 区 气温 是 影响 区 域内 植被 长 势 的 主导 因素 。 从 图 13b 
可 以 看 出 ,农业 生态 区 与 气温 呈正 相关 的 区 域 主要 分 布 在 河北 平原 .河南 东南 部 .安徽 北部 ,呈现 负 相 关 的 区 
域 主要 集中 于 山东 半岛 以 及 山东 省 西北 部 地 区 , 且 相 对 集中 。 华 北 农 业主 要 以 种 植 冬小麦 玉米、 棉花 为 主 ， 
冬小麦 生长 周期 长 ,上 且 在 主要 生长 期 多 抽取 地 下 水 灌溉 ,导致 区 域 植被 覆盖 度 与 降水 量 相 关 性 较 弱 ,甚至 呈 负 


http://www.ecologica.cn 


201702.00115v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


588 生 态 学 报 37 卷 


相关 关系 。 和 森林 生态 区 的 植被 覆盖 度 与 年 均 温和 年 降水 量 的 平均 相关 系数 为 0.03 和 -0.001 ,与 之 对 应 的 呈 
正 相关 区 域 面积 占 比 为 53% 和 49% ,表明 气温 是 影响 森林 生态 区 植被 覆盖 度 的 主导 因 了 于 。 和 森林 生态 区 植被 
履 盖 度 与 降水 呈现 负 相关 ,原因 在 于 降雨 意味 着 辐射 碱 少 和 温度 的 降低 ,从 而 影响 植被 生长 ”1 。 森 林 生态 
区 植被 覆盖 度 与 气温 因子 的 相关 性 有 区 域 差异 。 北 部 太行 山 ,燕山 山脉 的 植被 生长 更 多 受到 低温 的 影响 , 因 
此 与 气温 呈正 相关 的 区 域 面积 更 大 , 约 为 55% 。 山 东 半 岛 丘 陵 地 区 的 森林 植被 覆盖 度 与 气温 呈现 正 相 关 的 
区 域 占 比 为 35%。 山 东 半 岛 丘 陵 地 区 相 较 于 太行 山 .燕山 山脉 更 靠 南 ,植被 生长 所 需 的 热量 资源 更 丰富 。 热 
量 资源 较为 充足 的 地 区 , 增 温 对 植被 的 影响 会 变 弱 :” 。 太 行 山脉 地 区 的 降水 量 与 植被 覆盖 度 相 关 性 高 于 山 
东 半 岛 丘 陵 地 区 的 相关 性 ,燕山 地 区 的 相关 性 最 低 。 除 了 各 地 区 的 气候 地 形 .土壤 类 型 不 同 以 外 ,太行 山 山 
脉 的 植被 类 型 以 落叶 灌 和 草包 为 主 ,山东 半岛 丘陵 地 区 植被 以 温带 针 叶 林 和 阀 叶 林 为 主 ,燕山 山脉 的 植被 类 
型 以 灌 从 和 温带 阔 叶 林 为 主 。 不 同 的 植被 类 型 ,对 降水 的 响应 机 制 不 同 ,增强 了 相关 系数 的 空间 异 质 性 。 草 
原生 态 区 植被 覆盖 度 与 年 均 温和 年 降水 量 的 平均 相关 系数 为 -0.015 和 0.18 ,与 之 对 应 的 呈正 相关 的 区 域 面 
积分 别 占 分 区 面积 的 39% 和 92%. 说 明 降 水 量 是 影响 该 区 域 植 被 生长 的 主要 气候 因子 。 该 结论 与 张 戈 丽 和 称 
少 杰 ”的 结论 一 致 。 该 区 域内 的 草原 生态 区 位 于 内 蒙古 高 原 , 地 处 我 国 水 资源 相对 紧缺 地 区 ,因此 降水 对 
该 地 区 的 影响 要 高 于 气温 的 影响 。 
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图 13 植被 覆盖 度 与 降水 量 (a) 和 温度 (b) 的 相关 系数 
Fig.13 ”Coefficient of vegetation and temperature (a) and precipitation (b) 


3.4 气候 对 植被 覆盖 度 的 滞后 性 研究 

以 1982 一 2006 年 植被 覆盖 度 和 气象 数据 为 基础 ,分 别 计 算 生 长 季 3 月 一 11 月 的 月 植被 覆盖 度 和 月 降水 
量 月 均 温 的 偏 相关 系数 ,以 探讨 月 植被 覆盖 度 与 前 一 月 以 及 当月 气候 因子 的 相关 性 ,进而 得 出 植被 覆盖 度 对 
气候 因子 的 滞后 性 变化 。 

华北 平原 当月 植被 覆盖 度 与 当月 的 月 均 温 和 月 降水 量 的 平均 相关 系数 为 0.14 和 0.13( 图 14) ,与 之 对 应 
的 呈正 相关 的 区 域 面积 分 别 占 总 面积 的 88% 和 94% ,其 中 通过 P<0.05 检验 的 面积 占 比分 别 为 27.45% 和 
31.15%。 华 北平 原 月 植被 覆盖 度 与 前 一 月 的 月 均 温 和 月 降水 量 的 平均 相关 系数 为 0.008 和 0.227( 图 15), 与 
之 对 应 的 呈正 相关 的 区 域 面积 分 别 占 总 面积 的 47% 和 85%, 通 过 P<0.05 检验 的 面积 占 比分 别 为 8.01% 和 
54.65%。 总 体 而 言 ,植被 覆盖 度 与 前 一 个 月 的 月 均 温 的 相关 性 低 于 当月 均 温 的 相关 性 ,与 前 一 个 月 的 月 降水 
量 的 相关 性 高 于 当月 降水 量 的 相关 性 。 说 明 在 月 尺度 上 ,降水 量 对 整体 植被 覆盖 度 有 时 滞 效 应 。 植 被 覆盖 度 
与 前 一 月 气象 因子 关系 的 空间 分 布 具有 显著 纬度 地 带 性 :月 均 温 与 植被 覆盖 度 的 相关 系数 呈现 随 纬度 升 高 而 
降低 的 总 体 趋 势 ,月 降水 量 与 植被 覆盖 度 的 相关 系数 呈现 随 纬度 增高 而 升 高 的 趋势 。 总 体 月 尺度 的 相关 性 要 
远 高 于 年 际 水 平 。 

农业 生态 区 植被 覆盖 度 与 前 一 月 降水 量 的 平均 相关 系数 为 0.3 ,高 于 其 与 当月 降水 量 的 相关 系数 0.131。 
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其 与 前 一 月 的 月 均 温 的 平均 相关 性 为 -0.08, 低 于 其 与 当月 均 温 的 相关 性 0.109 ,这 表明 ,在 月 水 平 上 农业 生态 
区 植被 生长 对 降水 量 的 响应 存在 时 滞 效 应 ,与 月 均 温 则 没有 明显 的 时 滞 效 应 。 气 温 随 季节 变化 是 一 种 自然 规 
律 ,在 植物 生长 季 ,气温 是 植被 某 些 生长 期 的 关键 诱因 ,植被 遵循 气温 的 变化 进入 不 同 的 生长 期 ,因而 在 月 尺 
度 上 ,气温 对 植物 生长 不 存在 时 滞 效 应 。 森 林 生 态 区 植被 覆盖 度 与 前 一 月 的 月 均 温 的 平均 相关 系数 为 0.13 ， 
低 于 其 与 当月 的 数值 0.205 ,其 与 前 一 月 降水 量 的 相关 系数 为 0.12, 与 当月 降水 量 的 数值 0.11 基本 相同 ,整体 
而 言 不 存在 时 滞 效 应 。 这 与 穆 少 杰 :5 的 研究 结果 一 致 。 其 中 当月 的 月 均 温 与 月 植被 覆盖 度 的 相关 系数 呈正 
相关 的 面积 占 分 区 面积 的 93% ,范围 分 布 在 0 一 0.35 之 间 。 草 原生 态 区 植被 覆盖 度 与 前 一 月 的 月 均 温 的 平均 
相关 系数 为 0.24, 高 于 其 与 当月 平均 相关 性 0.046 ,其 与 前 一 月 降水 量 的 相关 系数 为 0.07, 低 于 其 与 当月 平均 
相关 性 系数 0.26。 这 表明 ,草原 生态 分 区 对 气温 的 响应 具有 时 滞 效 应 。 在 月 水 平 上 ,相关 性 结果 表现 出 一 定 
的 经 度 地 带 性 :从 西 到 东 ,植被 生长 与 降雨 的 相关 性 递减 ,同时 与 温度 的 相关 性 则 表现 出 递增 趋势 。 
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Fig.14 The coefficient of vegetation and temperature (a) and precipitation (b) in month scale 
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图 15 植被 覆盖 度 与 前 一 月 温度 (a) 和 降水 量 (b) 的 相关 系数 


Fig.1S The coefficient of vegetation and temperature (a) and precipitation (b) of the month bfore 


3.5 人 类 活动 对 植被 覆盖 度 影响 力 的 趋势 分 析 

为 探讨 人 类 活动 对 植被 纱 盖 度 的 影响 ,本文 用 残 差分 析 计 算 人 类 活动 对 植被 覆盖 度 的 影响 ,并 在 此 基础 
上 计算 人 类 活动 影响 的 变化 趋势 。 图 16 为 影响 植被 覆盖 度 变 化 的 人 类 活动 影响 力 的 变化 趋势 空间 分 布 图 。 
统计 结果 显示 ,人 类 活动 影响 力 的 整体 平均 变化 趋势 为 0.004/10a, 说 明 人 类 活动 影响 力 呈 现 继续 增强 的 趋 
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势 。 图 中 正 值 表示 人 类 活动 影响 力 继续 增强 , 负 值 表示 
影响 力 减 弱 。 从 图 中 可 见 , 河 北 北部 和 河南 中 东部 地 区 
人 类 活动 影响 力 在 持续 下 降 ,其 中 京 津 唐 地 区 减弱 速度 
最 快 。 以 北京 市 为 例 , 不 断 涌 人 的 人 口 加 快 了 城市 化 进 
程 ,中 心 城区 的 不 透水 面 建设 接近 饱和 ,导致 人 类 对 杆 
被 覆盖 度 变 化 的 影响 力 的 变 弱 。 山 东 省 西南 部 .安徽 和 
江苏 北部 地 区 人 类 活动 影响 力 呈 现 持续 增强 的 态势 。 
从 生态 学 角度 出 发 ,草原 生态 区 的 人 类 影响 力 增加 速率 
最 快 (0.005/10a), 农 业 生 态 区 的 影响 力 增加 最 慢 
(0.003/10a) 。 农 业 生 态 区 中 ,果树 和 粮食 混 种 区 的 人 
类 影响 增加 的 趋势 最 快 (0.005/10a) ,粮食 和 经 济 农 作 
物种 植 区 次 之 (0.004/10a) ,果树 种 植 区 最 慢 (0.003/ 115°00'00" 120°00'00"E 
10a) 。 同 时 ,结合 植被 覆盖 度 变化 趋势 图 ,统计 研究 区 ER 
之。 ”内 人 类 活动 的 整体 作用 力 为 -<5%( 其 中 ,植被 覆盖 度 呈 
区 ) 增加 趋势 的 区 域 的 相对 作用 力 为 正 值 ,植被 覆盖 度 变化 
趋势 呈 退 化 趋势 区 域 的 相对 作用 力 为 负 值 ,两 者 之 差 即 为 整体 相对 作用 力 ) ,说 明 人 类 活动 将 继续 以 负面 影 
响 的 作用 力 为 主 。 气 候 突变 前 植被 改善 区 的 人 类 活动 作用 力 低 于 植被 退化 区 的 相对 作用 力 ; 气 候 突 变 后 植被 
改善 区 的 人 类 活动 作用 力 高 于 退化 区 。 农 业 生 态 区 的 作用 力 一 直 呈 现 负面 作用 ,改善 区 的 相对 作用 力 减 弱 了 
4%。 草 原生 态 区 的 作用 力 呈 现 改 善 力 增强 的 趋势 ,幅度 达到 12%。 
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Fig.16 Dynamic of vegetation coverage 


二 ”4 结论 与 讨论 


4.1 结论 

本 文 基于 CIMMS-NDVI 和 MODIS-NDVI 遥感 数据 ,结合 1981 一 2013 年 气象 数据 ,利用 残 差分 析 .趋势 分 
局 析 、 相 关 性 分 析 ,揭示 了 华北 平原 近 30 年 气候 和 植被 覆盖 度 的 空间 变化 特征 ,并 对 气候 变化 对 各 生态 区 内 植 
人 被 覆盖 度 的 影响 进行 了 分 析 , 结 论 如 下 : 
rE (1)1981 一 2013 年 间 , 华 北平 原 气温 整体 呈现 显著 上 升 趋势 (0.20%C/10a,P<0.01) ,一 方面 全 球 气候 变 暖 
轧 ”背景 下 大 部 分 站 点 气温 整体 升 高 , 另 一 方面 春季 气温 的 升 高 .>15%C 等 温 线 控制 范围 的 扩大 和 年 均 温 0%C 等 值 
线 在 华北 平原 的 消失 ,也 是 区 域 平均 气温 升 高 的 诱因 。 其 中 ,15C 等 值 线 的 年 代 际 间 南 北 波动 范围 达到 2 个 
纬度 , 有 大 于 200km 的 北 移 。 华 北平 原 降 水 整体 呈现 显著 减少 趋势 (-1.75mm/10a,P>0.05) ,其 中 秋季 降水 
量 减少 过 快 ,加 上 400 一 600mm 降水 等 值 线 控制 范围 的 扩大 ,压缩 了 600 一 800mm 和 800 一 1000mm 降水 等 值 
线 的 范围 ,这 些 因素 都 造成 了 区 域 降水 量 的 减少 。 四 季 气 候 倾向 率 的 特征 变化 敏感 区 域 主要 位 于 35° 一 39°N 
之 间 。1991 一 1994 年 为 华北 平原 气候 的 突变 时 期 。 

(2) 华北 平原 植被 覆盖 度 空间 分 布 具有 显著 的 地 区 差异 性 ,整体 上 呈现 上 升 趋势 ,呈现 增加 趋势 的 面积 
占 比 为 55% 。 草 原生 态 区 的 年 均 植被 覆盖 度 较 低 ,森林 生态 区 和 农业 生态 区 相对 较 高 。 从 植被 覆盖 度 和 植 
被 覆盖 度 变 化 趋势 的 空间 分 布 来 看 ,人 类 活动 不 仅 加 速 了 区 域 植被 覆盖 度 的 降低 ,也 加 剧 了 降低 速率 的 变 快 。 
在 基于 年 代 际 的 时 间 尺 度 研究 中 ,气温 的 增加 ,对 华北 平原 农业 生态 区 的 影响 体现 为 植被 覆盖 度 的 稳定 增加 ， 
但 是 速率 在 不 断 降低 ;草原 生态 区 植被 覆盖 度 的 变化 趋势 在 一 定 程度 上 与 气温 的 变化 趋势 一 直 。 当 气温 偏 高 
降雨 偏 多 时 ;森林 生态 区 的 植被 覆盖 度 呈 减少 趋势 ,其 余 的 水 热 组 合 方式 植被 覆盖 度 均 呈 现 增加 趋势 。 

(3) 总 体 上 ,研究 区 月 尺度 植被 覆盖 度 与 气候 因子 的 相关 性 高 于 年 尺度 上 的 值 。 植 被 覆盖 度 与 年 降水 量 
的 偏 相 关 性 高 于 其 与 年 均 温 的 偏 相关 性 ,但 是 存在 明显 的 区 域 差 异 。 年 均 温 对 农业 生态 区 和 和 森林 生态 区 的 植 
被 覆盖 度 的 影响 更 大 ,草原 生态 区 的 植被 覆盖 度 对 年 降水 量 的 依赖 性 更 大 。 在 月 尺度 上 年 降水 量 对 植被 的 影 
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响 具 有 时 滞 效 应 。 气 温 对 草原 生态 区 具有 时 滞 效 应 ,降雨 对 农业 生态 区 具有 时 清 效 应 。 

(4) 华 北平 原 干 热 化 的 气候 突变 降低 了 植被 覆盖 度 的 增加 趋势 。 在 气候 突变 前 ,增加 速率 最 快 的 是 农田 
生态 区 ;气候 突变 后 草原 生态 区 植被 覆盖 度 增 加 速率 最 快 。 从 植被 覆盖 度 高 低 的 角度 出 发 ,气候 突变 对 草原 
生态 区 影响 最 大 ;从 变化 速率 角度 出 发 ,对 农田 生态 区 的 影响 最 大 。 虽 然 气候 突变 后 ,人 类 活动 在 植被 覆盖 度 
的 改善 区 的 相对 作用 力 高 于 退化 区 ,但 是 就 整体 而 言 , 人 类 活动 的 影响 力 还 在 持续 增强 ,上 且 呈 现 出 在 退化 区 的 
作用 力 高 于 改善 区 的 趋势 。 

4.2 讨论 

本 文 从 生态 学 角度 出 发 ,对 不 同 生态 区 植被 覆盖 度 的 分 布 与 变化 趋势 进行 了 分 析 讨 论 ,总 结 了 基于 生态 
学 分 区 的 植被 覆盖 的 分 布 和 变化 规律 ,探讨 了 植被 覆盖 度 与 气候 因子 的 关系 。 相 较 于 传统 统计 分 析 结 果 , 相 
关系 数 有 一 些 不 同 。 与 把 整个 研究 区 域 作 为 统计 单元 的 研究 方法 相 比 ,基于 像 元 的 小 样本 统计 能 更 好 的 体现 
NDVI 和 和 气候 因子 的 空间 差异 性 。 在 统计 年 限 有 限 的 情况 下 ,显著 的 空间 差异 可 能 会 影响 小 样本 的 统计 分 析 
结果 ,掩饰 了 气候 因子 与 NDVI 的 相关 性 。 另 外 ,植被 生长 季 的 有 效 降 水 跟 NDVI 的 相关 性 可 能 更 高 ,此 次 研 
究 没 有 排除 非 生 长 季 的 无 关 水 影响 。 

水 热 组 合 是 决定 植被 生长 的 重要 非 生物 因素 。 虽 然 本 研究 利用 不 同年 代 际 的 水 热 组 合 对 植被 覆盖 度 的 
响应 进行 了 粗略 探讨 ,但 是 水 热 组 合 的 形式 多 样 , 且 不 同 植物 种 类 对 水 热 组 合 的 胁迫 响应 机 制 也 各 不 相同 。 
适宜 的 水 热 组 合 加 速 植物 的 生长 ,不 适宜 的 水 热 组 合 将 对 植被 生长 产生 抑制 作用 。 因 此 ,在 全 球 气候 不 断 变 
暧 ,厄尔尼诺 现象 频 发 的 大 背景 下 ,研究 植被 对 不 同 水 热 组 合 的 响应 机 理 ,探讨 不 同 植被 类 型 气候 变化 的 响应 
和 分 布 , 可 以 做 为 下 一 步 研究 内 容 。 

除了 气候 因子 ,人 类 活动 也 是 影响 植被 覆盖 度 变 化 的 重要 驱动 因素 ,虽然 本 文 对 人 类 活动 的 整体 影响 做 
了 简单 分 析 ,但 是 人 类 活动 对 植被 覆盖 度 的 影响 机 理 错综复杂 。 特 别 是 进入 21 世纪 ,人 类 活动 强度 和 范围 逐 
渐 扩 大 ,包括 人 口 增 长 城镇 化 发 展 能 源 消耗 结构 调整 .退耕 还 林 .退耕 还 牧 .农业 结构 的 调整 等 ,都 会 对 区 域 
植被 覆盖 度 产 生 严 重 影响 。 因 此 ,如 何 细 化 人 类 活动 对 植被 覆盖 度 的 影响 ,定位 重要 的 人 类 活动 因子 ,消除 影 
响 植 被 变化 的 不 利 因素 ,将 对 植被 覆盖 度 的 变化 产生 积极 影响 ,对 于 区 域 生态 安全 也 有 积极 意义 。 
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